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Abstract 

The behaviour of the magnetic field inside the superconducting quark mat- 
ter core of a neutron star is investigated in the framework of the Ginzburg- 
Landau theory. We take into account the simultaneous coupling of the diquark 
condensate field to the usual magnetic and to the gluomagnetic gauge fields. 
We solve the problem for three different physical situations: a semi-infinite re- 
gion with a planar boundary, a spherical region, and a cylindrical region. We 
show that Meissner currents near the quark core boundary effectively screen 
the external static magnetic field. 
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1 Introduction 



Recently, possible formation of diquark condensates in QCD at finite density has 
been re- investigated in series of papers following Refs. [£[]], 0. It has been shown 
that in chiral quark models with non-perturbative 4-point interaction motivated 
from instantons || or non-perturbative gluon propagators ||, || the anomalous 
quark pair amplitudes in the color antitriplet channel can be very large: of the 
order of 100 MeV. Therefore, one expects the diquark condensate to dominate the 
physics at densities beyond the deconfmement /chiral restoration transition density 
and below the critical temperature (of the order of 50 MeV). Various phases are 
possible. The so called two-flavor (2SC) and three-flavor (3SC) phases allow for 
unpaired quarks of one color. It has been also found ||, 0] that there can exist a 
color-flavor locked (CFL) phase for not too large strange quark masses ||, where 
color superconductivity is complete in the sense that diquark condensation results 
in a pairing gap for the quarks of all three flavors and colors. 

The high-density phases of QCD at low temperatures are relevant for the ex- 
planation of phenomena in rotating massive compact stars which might manifest 
themselves as pulsars. Physical properties of these objects (once being measured) 
could constrain our hypotheses about the state of matter at the extremes of den- 
sities. In contrast to the situation for the cooling behaviour of compact stars 0, 
[10(1, where the CFL phase is dramatically different from the 2SC and 3SC phases, 



we don't expect qualitative changes of the magnetic field structure for these phases. 
Therefore, below we shall restrict ourselves to the discussion of the simpler two- 
flavor theory first. Bailin and Love |TT| used a perturbative gluon propagator which 
yielded a very small pairing gap and they concluded that quark matter is a type I 
superconductor, which expells the magnetic field of a neutron star within time-scales 
of 10 4 years. If their arguments would hold in general, the observation of life-times 
for magnetic fields as large as 10 7 years |T2]], |TB| would exclude the occurence of 



an extended superconducting quark matter core in pulsars. These estimates are not 
valid for the case of diquark condensates characterized by large quark gaps. Besides, 
in Ref. |T3[] the authors found that within recent non-perturbative approaches for 



the effective quark interaction that allow for large pairing gaps the quark condensate 
forms a type II superconductor. Consequently for the magnetic field H < H c \ there 
exists a Meissner effect and for H C 2 > H > H c \ the magnetic field can penetrate 
into the quark core in quantized flux tubes. However, they have not considered 
in that paper the simultaneous coupling of the quark fields to the magnetic and 
gluomagnetic gauge fields. 

Though color and ordinary electromagnetism are broken in a color superconduc- 
tor, there is a linear combination of the photon and the gluon that remains massless. 
The authors of Ref. (15| have considered the problem of the presence of magnetic 
fields inside color superconducting quark matter taking into account the possibility 
of the so called "rotated electromagnetism". They came to the conclusion that there 
is no Meissner effect and the external static homogeneous magnetic field can pene- 
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trate into superconducting quark matter because in their case it obeys the sourceless 
Maxwell equations. To our opinion this result is obtained when one does not pose 
correct boundary conditions for the fields. Obviously it is energetically favorable to 
expell the magnetic field rather than to allow its penetration inside the supercon- 
ducting matter. Using for the description of the diquark condensate interacting with 
two gauge fields the same model as in Refs. |J, f| and we will show that the 



presence of the massless excitation in the spectrum does not prevent the Meissner 
currents to effectively screen the static external magnetic field. 



In Ref. [IB | two of us have derived the Ginzburg-Landau equations of motion 



for the diquark condensate placed in static magnetic and gluomagnetic fields 

e -> q 

■o A + 4= 
3 y/3 



ad p + (3{d p d*)d p + 7 (i y -\1 +4r G 8 ) 2 d p = 0, (1) 



where d p is the order parameter, a = t dn/dE, (3 = (dn/dE) 7 £(3)(7rT c )~ 2 /8, 7 = p 2 F 
/3/(6yU 2 ), dn/dE = p F fi/n 2 , t = (T — T c )/T c , T c being the critical temperature, p F - 
the quark Fermi momentum, and for the gauge fields 



sma(dpVd p - d p Vd p ) 
2qd 2 



A ? rot rot A + sin a A— % 7T~n — 1" sm a cos a Gs, (2) 



2 -* 2 -» cos a(d p V d* — d*V d p ) -» 
A„rot rot G 8 + cos a G8= — « ^r - ^ — 1" sin a cos a A • (3) 

These equations introduce a "new" charge of the diquark pair q = ^/rfe 2 + g 2 Ps, 
Ps — 1/ \/3, and for the diquark condensate with paired blue-green and green-blue ud 
quarks one has r] = 1/ y/3. The penetration depth of the magnetic and gluomagnetic 
fields X q and the mixing angle a are given by 

K 1 = Q d yfir> cosa = ^w + g* ■ ^ 

At neutron star densities gluons are strongly coupled (g 2 /An ~ 1) whereas photons 
are weakly coupled (e 2 /47r = 1/137), so that a ~ rje/g is small. For g 2 /At: ~ 1 we 
get a ~ 1/20. The new charge q is by an order of magnitude larger than e/\/3. 
Please notice also that since red quarks are normal in the 2SC and 3SC phases, 



there exist the corresponding normal currents j^(A) = —Tl^A" and j^(Gg) = 
— Il^Gg which however do not contribute in the static limit under consideration 



to the above Ginzburg - Landau equations, cf. |Tj . Thus, the qualitative behavior 
of the static magnetic field for all three 2SC, 3SC, and CFL phases is the same. 

We will solve the Ginzburg-Landau equations ([[]), (fj), ([|) for the case of a 
homogeneous external magnetic field for three types of superconducting regions: 
a) a semi-infinite region with a planar boundary, b) a cylindrical region and c) a 
spherical region. The ) and c) simulate the behavior of the magnetic field in 

quark cores of massive neutron stars. A discussion of all the three cases might be 
interesting in connection with the expectation that the slab, the rod and the droplet 
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structures may exist within the mixed quark-hadron phase of the neutron stars, cf. 

hi, irran. 



We assume a sharp boundary between the quark and hadron matter since the 
diffusion boundary layer is thin, of the order of the size of the confinement radius 



(~ 0.2 0.4 fm), and we suppose that the coherence length = J"y/(—2a) is not 
less than this value and the magnetic and gluomagnetic field penetration depth X q 
is somewhat larger than the confinement radius. Also we assume that the size of 
the quark region is much larger than all mentioned lengths. 

2 Solution of Ginzburg-Landau equations 

Let us rewrite equations (Q) and (|3]) for a homogeneous superconducting matter 
region (being either a type I superconductor, or a type II superconductor for H < 
H c i) in the following form 

A 2 rot rot A + sin 2 a A— sin a cos a Gs, (5) 

A 2 rot rot Gs + cos 2 a Gs= sin a cos a A ■ (6) 
The field Gs is defined from (H) as follows 

-> A 2 , rot rot A + sin 2 a A 

G 8 = : • 7 

sin a cos a 

From equations (|5]) and (0) we obtain the relation 

rot rot (^8= — cot a rot rot A ■ (8) 
Substitution of Gs from ([?]) into (||) yields 

rot rot (A 2 rot rot A + A) = 0. (9) 
Introducing the new function M, 

M= rotrot A , (10) 

we obtain 

A 2 , rotrot M + M = 0. (11) 

Thus the vector potential A can be determined by simultaneous solution of equations 
( |10"D and (|TT|), whereas the gluonic potential Gs is found from (0). 



For the solution of equations (^), (|T0|) and (|TTD we also need appropriate bound- 
ary conditions. At the quark-hadronic matter boundary we require the continuity 
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of the magnetic field and the vanishing of the gluon potential (Gs— 0) due to gluon 
confinement. Also the potential Gs and the magnetic induction can't be infinite 
within the region of their existence. As we shall see below, these conditions are 
sufficient for a unique determination of the magnetic and gluomagnetic fields inside 
the quark matter. 

Equations (|TJ), (|2|) and (^) have an obvious solution for a homogeneous (Vrf p = 0) 
and infinite superconductor (A = 0, Gs= 0) 

d 2 = -a/ (3 > 0. (12) 

This solution motivates the possibility of existence of the complete Meissner effect 

for both the fields A and Gs inside the quark superconductor. It corresponds to 
the absolute minimum value of the free energy f = f n — a 2 / (2(3) , where f n is 
the free energy of the normal quark matter fllTf , ||14|| . The presence of the fields 

A 7^ 0, Gsi^ inside the quark region would increase the free energy. 



As we have mentioned, the estimates [14], [16] have demonstrated that color 
superconductors are type II superconductors. Indeed, the Ginzburg - Landau pa- 
rameter k = \ q /l^ = Vfi/do) > 3. For type II superconductors one can drop the 

fields A and Gs in the Ginzburg - Landau equation ([!]) arriving at the solution (|12"D , 
since the coherence length 1% is smaller than the penetration depth of the magnetic 
and gluomagnetic fields X q . Then in equations of motion (fj) and (||) the penetration 
depth of the magnetic and gluomagnetic fields \ q can be put constant. This simpli- 
fies solution of the Ginzburg - Landau equations very much. For further analytical 
treatment of the problem we will assume that we deal with a type II superconductor 
although our main conclusion on the existence of the Meissner effect is quite general. 

2.1 Planar boundary 

Let us first consider a semi-infinite region of superconducting quark matter for x < 
with a planar boundary, which coincides with the zy plane. The external static ho- 
mogeneous magnetic field H is directed along the z-axis, the external vector potential 
A is aligned in the y direction A y (x) = Hx. Then the internal potentials A and Gs 
also have only y components: A y = A y (x) and Gs y = Gs y (x), div A = 0, div Gs= 0. 



Simultaneous solution of equations fllCf ) and (|Tl"|) determines the electromagnetic 
vector potential A y as follows 

x 

A y = ci exp(— ) + c 2 x + c 3 . (13) 
\ 

We put c 2 = because otherwise A y — > oo and G 8y — > oo for x — » oo that would 
lead to a complete destruction of the condensate. We search for an energetically 
favorable unique solution of the problem satisfying the above mentioned boundary 
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conditions. Thus we further put 

A y = ciexp(^) + c 3 . (14) 

Substitution of the solution ( |14|) into equation (0) yields the gluonic potential Gs y , 

x 

Gs y = — cot a c\ exp( — ) + tan a C3 . (15) 

\ 

The boundary condition G$ y (x = 0) = determines the constant C3 = cicot 2 a . 
For the vector potential of the magnetic field we obtain 

A y = ci(exp(— ) + cot 2 a). (16) 
\ 

We can determine c\ from the remaining boundary condition dA y (x = 0)/dx = 
which yields c\ = HX q . Finally the potentials render 

A y = H\ q (exp(^) + cot 2 a) , (17) 

x 

Gs y = — cot a H Aq(exp(— ) — 1) . (18) 

\ 

Then, for the magnetic induction B z = dA y /dx and for the gluonic field K z = 
dG 8y / dx inside the quark superconductor we have the following expressions 

B=Hexp(^) e z , (19) 
K= - cot a B l z . (20) 



Ref. [15[ introduced the "rotated" fields B x and B y , 

B x =-sina B +cosa K, (21) 

B a = cos a B+ sin a K . (22) 
It is easy to see that solutions (|T9| ) and (|20|) yield B y z = and B x z = —(if/sin a) exp(x/ \ q ) 

Therefore there is no B y field inside the superconducting quark matter and the B x z 
field is expelled from the color superconductor. Consequently, in contradiction with 



the statement of Ref. [15| there is a Meissner effect for the color superconductors. 
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2.2 Cylindrical structures 

Now we shall consider a cylindrical region of superconducting quark matter of radius 
a , whose axis coincides with z axis of cylindrical coordinates. Such a situation may 
occur in the mixed phase where the rods of quark matter imbedded in the hadron 
matter are possible configurations. The external homogeneous magnetic field H 
is directed along the z axis, the external vector potential A is aligned in the (p - 
direction = Hr/2 . The internal vector potentials A and Gs have only tp - 
components A v {r) and G^{r) . Thus, equation ([LI]) acquires the form 



d 2 M in 1 dM tr 



dr 2 



+ - 



r dr 



. 



(23) 



The solution of this equation is M 9 = — Cil\{r/\ q ), where Ii(r/X q ) is the corre- 
sponding modified Bessel function. Consequently, equation (|i~0"D can be written in 
the following form 



d l A„ 



+ 



1 dA u 



dr 2 ' r dr 
For the vector potential A v we obtain 



A 

r' 



(24) 



Mr) 



T 

ciXpi(-) + c 2 r. 



(25) 



Then for the potential from (|^) we find the following expression 



G 8ip {r) 



cot a C\ Ag 1\ ( — ) + tan a c 2 r. 



(26) 



The magnetic induction B z inside the superconductor is given by 



B z {r) = ciA 9 / ( — ) + 2c 2 . 

We determine c 2 from the boundary condition G 8lfi (a) ■ 
Consequently, final expressions for the potentials are 



and Ci from BAa) 



(27) 
= H. 



Mr) 



P(a/X q ,a) 



h{— ) + cot 2 a 

A n a A 



HX q cot a 



P(a/X g ,a) [ X g 
and for the corresponding fields we get 



Mf ) - -Mf ) 

a Aq 



(28) 
(29) 



H 



Bz r P{a/X q ,a) [' ux X q 



J (^) + 2cot 2 a^/ 1 (-^; 



a 



X 



i J 



(30) 
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KJr) 



H cot a 
P(a/X q ,a) 



a 



A 



where P(a/X q ,a) is given by 



P(a/X q ,a) 



X n 



+ 2^cot 2 a/ 1 (-^) . 
a A„ 



Thus, we obtain the following final expressions for the rotated fields 



B v , 



B\ 



2X q H cot a 
aP(a/X q , a) sin a; 



P(a/A, 



sin a 




(31) 
(32) 

(33) 
(34) 



We notice that the field B y is homogeneous inside the quark matter. Thus one may 
expect that in presence of the magnetic field the slab structures are energetically 
favorable starting with a smaller quark fraction volume than in absence of the mag- 
netic field since, as we have argued, the magnetic field is expelled from the slabs 
and it penetrates the cylindres. In application to the description of the quark core 
of the neutron star, the discussion of the cylindrical case has no meaning and one 
should further consider the case of a spherical geometry. 



2.3 Spherical case 

We now assume that the neutron star possesses a spherical core of radius a consisting 
of color superconducting quark matter. The applied homogeneous magnetic field H 

is directed along the z axis. The functions M, A and Gg have only tp - components: 
M v (r, A v (r,$) and G 8ip (r,{}). For the solution of the equation flllD we use the 
ansatz M v (r, i?) = f(r) sin$. Then the equation ( |TTD can be written in spherical 
coordinates as 

d 2 f 2df (2 1\ 



dr 2 r dr 

The solution of equation (]35f) which tends to zero at the centre of the core renders 

f(r) = —Ay ), (36) 

' A q 



where 



J(-) = sinh(-)--cosh(-). (37) 

A q A q A q A q 



For the solution of equation (1C) we use the ansatz Ay,(r, i?) = g(r) siw&. Then 
equation (|1^ ) acquires the form 

d 2 q 2 do 2q D , r , . . 

dr A r dr r z r z X„ 
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The general solution of equation (|H]) in spherical coordinates is given by 



DM 



J(— ) sin$ + C\r sin-i?. 
A n 



(39) 



We also find the general solution for the gluonic potential from (|7]) 

D\ 2 r 

GgJr, $) = — cot a — tt'J(t-) sint? + tana C\T sin •&. 
r z X„ 



(40) 



The constant c\ is determined from the boundary condition G 8lfi (a, i?) = as q = 
-DA 2 cot 2 aJ(a/X q )/a 3 . Then expressions for the vector potential and for the 
gluonic potential are written as 



j ) + cot 2 a -J - 
Ag a d A g 



cot a- 



DA 



2 r 



A 9 a d A 9 j 



sin??, 
sin0. 



(41) 



Therefore, the expressions for the radial component B ir and the tangential compo- 
nent acquire the form 



where 



Ji( 



2DA 



a , r , 2 w a 



-J(-) + cofa J(— ) 
r d A 9 A 9 



-DA 2 la 3 



r 



3 J 1 (-)-2cot 2 aJ(-^) 

6 A„ An 



cost?, 
sin0 , 



(42) 



A„ 



:i+ A| )smh( A: 



Aq A q 



(43) 

Accordingly, the expressions for the radial component K ir and the tangential com- 
ponent Ki$ of the internal gluomagnetic field render 



K ir {r,$) 



cot a- 



2DA 2 la 3 



— cot a- 



DXl 



a. 



COS0, 



(44) 



A„ 



sin 



The solutions of the Maxwell equations outside the sphere are 



2m 



., j COS0, 



(45) 



sin0. 
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We can find D and m from the boundary conditions B ir (a, i?) = B er (a, i?), B^(a, ■&) 
B e #(a,$). Thus, we obtain 



D = 



m 



3Hc 



2iVsinh(a/Ag)' 
Ha 3 



(46) 



2N 



1 + 3^-3^ coth^ 
a 2 a A„ 



where iV is given by 



^f-^coth^- 
cr a 



Xn 



N = 1 - 3 cot 2 a 

The final expressions for the internal magnetic fields have the form 

= - a^smhd/A,) [^ J( A; ) + COt aJ( A;r° S ^ 
3/7 A 2 r« 3 . 

Bj,,i — 



(47) 



(48) 



Ji( — ) - 2 cot 2 aJ(— ) 



A 



A n 



sini?. 



2a 2 iVsinh(a/A g ) 

The corresponding internal components of gluomagnetic field are given by 

3 if A 2 cot a 



K, 



3 

a r 



a 2 Nsinh(a/X q ) 

3HX 2 q cot a 
2a 2 Nsmh(a/X q) 



A 



A 



Q J 



Jl( +2 J - 



A n 



A, 



cost?, 
sini?. 



(49) 



We obtain therefore the following internal components of the rotated sourceless field 



Bl 



B y — 



3HX q cot a 
Na 2 sin a 
3HX 2 cot a 
Na 2 sin a 



a a 
— coth 1 

Xq Xq 

a a 

— coth 

L X q X q 



cos$, 



sin$. 



(50) 



It is to be noticed that the components of the rotated field B y depend only on the 
polar coordinate ■&. The internal components of the rotated massive field B x are 



3#A 2 a , r . 

Bf r = — — \. - r J — COST?, 

Afr 3 sinasinha/A, X q 



(51) 



BZ 



3HX 2 q a 



2Nr 3 smasmh(a/Xq) X 



Ji(— ) sini?. 
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In a neutron star, the radius of the superconducting quark core a is much larger 
than the penetration depth X q . In this limit the unbroken rotated field components 
take the form 



3H cot aA„ , 

B? = ^costf, (52) 

sin a a 

3H cot a \ n 



sin a a 

In the very same limit, for the components of the rotated massive field B x at the 
surface of the quark core we obtain 

B? r (a,0) = ^costf, (53) 

sin a a 

B&(a,0) = — sintf. 

sin a 

Therefore, the rotated massive field component at the surface of the quark core 
Bf# (a, 7?) is larger than by the factor a/\ q cot a and the B* r (a, •&) and Bf r com- 
ponents are of the same order of magnitude. We notice that there is Meissner effect 

for the rotated field B x . Let us now estimate the rotated field B y inside the star. 
For values of the external field strength H = 10 12 G, the penetration depth X q = 1.7 
fm, the mixing angle a = 0.05 and the radius of the quark core a = 1 km we obtain 

B y = 2- 10~ 4 G. Thus, we can safely neglect the rotated unbroken field B v . There- 
fore we again insist that there is a Meissner effect, and the applied external static 
magnetic field is thereby almost completely screened within the spherical geometry. 
Only a tiny fraction of the field can penetrate into the superconducting quark cores 
of the neutron stars. 



3 Conclusion 

We have investigated the behaviour of color superconducting quark matter in an 
external static homogeneous magnetic field and could show that color and electric 
Meissner currents exist. For this purpose we have solved the Ginzburg-Landau 
equations for three types of superconducting regions: a) a semi-infinite region with 
planar boundary, b) a cylindrical region and c) a spherical region. We have obtained 
analytic expressions for magnetic and gluomagnetic fields inside the quark matter 
for all three cases. In application to the quark cores of massive neutron stars we 

have shown that one can neglect the rotated field B v inside the superconducting 
quark core since the Meissner currents effectively screen the applied external static 



magnetic field. These results confirm the physical situation discussed in |T6j . 
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